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Abstract 

The astrophysical factor S pp (0) for the solar proton burning p + p — > D + e + 
+ v e is recalculated in the relativistic field theory model of the deuteron (RFMD). 
We obtain S , pp (0) = 4.08 x 10 -25 MeVb which agrees good with the recommended 
value 5pp(0) = 4.00 x 10 _25 MeVb. The amplitude of low-energy elastic proton- 
proton (pp) scattering in the x So-state with the Coulomb repulsion contributing 
to the amplitude of the solar proton burning is described in terms of the S-wave 
scattering length and the effective range. This takes away the problem pointed out 
by Bahcall and Kamionkowski (Nucl. Phys. A625 (1997) 893) that in the RFMD 
one cannot describe low-energy elastic pp scattering with the Coulomb repulsion 
in agreement with low-energy nuclear phenomenology. The cross section for the 
neutrino disintegration of the deuteron u e + D — ► e _ + p + p is calculated with 
respect to 5 PP (0) for neutrino energies up to E Uc < lOMeV. The results can be 
used for the analysis of the data which will be obtained in the experiments planned 
by SNO. The astrophysical factor S pep (0) for the process p + e~+p^i/ e + D (or 
pep-process) is calculated relative to <S' PP (0) in complete agreement with the result 
obtained by Bahcall and May (ApJ. 155 (1969) 501). 
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1 Introduction 

The weak nuclear process p + p^D + c + + z/ c , the solar proton burning or proton-proton 
(pp) fusion, plays an important role in Astrophysics [1,2]. It gives start for the p-p chain of 
nucleosynthesis in the Sun and the main-sequence stars [1,2]. In the Standard Solar Model 
(SSM) [3] the total (or bolometric) luminosity of the Sun L & = (3.846 ± 0.008) x 10 26 W 
is normalized to the astrophysical factor 5 , pp (0) for pp fusion. The recommended value 
5p P (0) = 4.00 x 10 -25 MeVb [4] has been found by averaging over the results obtained in 
the Potential model approach (PMA) [5,6] and the Effective Field Theory (EFT) approach 
[7,8]. However, as has been shown recently in Ref. [9] the inverse and forward helioseismic 
approach indicate the higher values of S pp (0) seem more favoured, for example, 5 , pp (0) = 
4.20 x 10~ 25 MeVb and higher [9]. Of course, accounting for the experimental errors the 
recommended value does not contradict the result obtained in Ref. [9]. 

In Refs. [10-13] we have developed a relativistic field theory model of the deuteron 
(RFMD). In turn, in Ref. [14] we have suggested a modified version of the RFMD which 
is not well defined due to a violation of Lorentz invariance of the effective four-nucleon 
interaction describing N + N — > N + N transitions. This violation has turned out to be 
incompatible with a dominance of one-nucleon loop anomalies which are Lorentz covari- 
ant. Thereby, the astrophysical factor S , pp (0) calculated in the modified version of the 
RFMD [14] and enhanced by a factor of 1.4 with respect to the recommended value [4] 
is not good established. This result demands the confirmation within the original RFMD 
[10-13] by using the technique expounded in Ref. [13]. 

As has been shown in Ref. [12] the RFMD is motivated by QCD. The deuteron ap- 
pears in the nuclear phase of QCD as a neutron-proton collective excitation - a Cooper 
np-pair induced by a phenomenological local four-nucleon interaction. Strong low-energy 
interactions of the deuteron coupled to itself and other particles are described in terms of 
one-nucleon loop exchanges. The one-nucleon loop exchanges allow to transfer nuclear 
flavours from an initial to a final nuclear state by a minimal way and to take into account 
contributions of nucleon-loop anomalies determined completely by one-nucleon loop di- 
agrams. The dominance of contributions of nucleon-loop anomalies has been justified in 
the large Nq expansion, where Nc is the number of quark colours [13]. Unlike the PMA 
and the EFT approach the RFMD takes into account non-perturbative contributions of 
high-energy (short-distance) fluctuations of virtual nucleon (JV) and anti-nucleon (N) 
fields, NN fluctuations, in the form of one-nucleon loop anomalies. In accord the analy- 
sis carried out in Refs. [15] nucleon-loop anomalies can be interpreted as non-perturbative 
contributions of the nucleon Dirac sea. The description of one-nucleon loop anomalies 
goes beyond the scope of both the PMA and the EFT approach due to the absence in 
these approaches anti-nucleon degrees of freedom related to the nucleon Dirac sea. How- 
ever, one should notice that in low-energy nuclear physics the nucleon Dirac sea cannot 
be ignored fully [16]. For example, high-energy NN fluctuations of the nucleon Dirac 
sea polarized by the nuclear medium decrease the scalar nuclear density in the nuclear 
interior of finite nuclei by 15% [16]. This effect has been obtained within quantum field 
theoretic approaches in terms of one-nucleon loop exchanges. 

In this paper we revise the value of <S , PP (0) obtained in Ref. [14]. For this aim we apply 
the technique developed in the RFMD [13] for the description of contributions of low- 
energy elastic nucleon-nucleon scattering in the 1 So-state to amplitudes of electromagnetic 



and weak nuclear processes. This technique implies the summation of an infinite series of 
one-nucleon loop diagrams and the evaluation of the result of the summation in leading 
order in large Nq expansion [13]. The application of this method to the evaluation of the 
cross sections for the anti-neutrino disintegration of the deuteron induced by charged i> e 
+ D^e + + n + n and neutral z/ c + D^z/ c + n + p weak currents gave the results 
agreeing good with the experimental data. The reaction v c + D — > e + + n + n is, in the 
sense of charge independence of weak interaction strength, equivalent to the reaction p + 
p — > D + e + + v e . Therefore, the application of the same technique to the description of 
the reaction p + p^D + e + + z/ c should give a result of a good confidence level. 

The paper is organized as follows. In Sect. 2 we evaluate the amplitude of the solar 
proton burning. We show that the contribution of low-energy elastic pp scattering in the 
1 So-state with the Coulomb repulsion is described in agreement with low-energy nuclear 
phenomenology in terms of the S-wave scattering length and the effective range. This 
takes away the problem pointed out by Bahcall and Kamionkowski [17] that in the RFMD 
one cannot describe low-energy elastic pp scattering with the Coulomb repulsion in agree- 
ment with low-energy nuclear phenomenology. In Sect. 3 we evaluate the astrophysical 
factor for the solar proton burning and obtain the value 5 , pp (0) = 4.08 x 10~ 25 MeVb 
agreeing good with the recommended one 5 , pp (0) = 4.00 x 10~ 25 MeVb. In Sect. 4 we 
evaluate the cross section for the neutrino disintegration of the deuteron v c + D — > e~ + 
p + p caused by the charged weak current with respect to S f pp (0). In Sect. 5 we adduce 
the evaluation of the astrophysical factor S , pep (0) of the reaction p + e~+p^D + z/ c 
or pep-process relative to S , pp (0). In the Conclusion we discuss the obtained results. 

2 Amplitude of solar proton burning and low— energy 
elastic proton— proton scattering 

For the description of low-energy transitions N + N— > N + N in the reactions n + p — > 
D + 7, 7 + D^n + p, z/ c + D— >e + + n + n and p + p — *■ D + e + + u c , where nucleons 
are in the 1 So-state, we apply the effective local four-nucleon interactions [11-13]: 

C^ NN (x) = G nNN {[n(x) lfll 5 p c (x)]lpc(x)Yl 5 n(x)] 

1 _ 1 , _ 

+ 2 \n(x)-i^n c (x)}\n c (x)^^n(x)} + - [p(x)-f^ 5 p c (x)}[p c (x)^ 5 p(x)} 

+ ( 7m7 5 ® 7 V^7 5 ®7 5 )}, (2.1) 

where n(x) and p(x) are the operators of the neutron and the proton interpolating fields, 
n c (x) = Cn (x) and so on, then C is a charge conjugation matrix and T is a transposition. 
The effective coupling constant G^nn is defined by [11-13] 

G.nn = f§| ~ ^ = 3-27 x 10- 3 MeV" 2 , (2.2) 

4M; M N 

where (JVnn = 13.4 is the coupling constant of the 7rNN interaction, M n = 135 MeV is 
the pion mass, M p = M n = Mn = 940 MeV is the mass of the proton and the neutron 
neglecting the electromagnetic mass difference, which is taken into account only for the 
calculation of the phase volumes of the final states of the reactions p + p^D + e + + 



u e , v e + D — > e~ + p + p and p + e~ + p — > D + u e , and a np = (—23.75 ± 0.01) fm is 
the S-wave scattering length of np scattering in the 1 S -state. 

The effective Lagrangian for the low-energy nuclear transition p + p— >D + e + + z/ e 
has been calculated in Ref. [12] and reads 

£pp^dc+,» = -ig A G„ NN M^^Dl(x) [p c (x) 7 5 p(x)] [^(^(l - 7 5 )V>e(z)]- (2.3) 

where Gy = Gp cos$c with Gp = 1.166 x 10~ n MeV~ 2 and $c are the Fermi weak 
coupling constant and the Cabibbo angle cos^c — 0.975, g A = 1.2670 ± 0.0035 [18] 
and gy is a phenomenological coupling constant of the RFMD related to the electric 
quadrupole moment of the deuteron Q D = 0.286 fm 2 [11]: g v = 2ir 2 Q D M^. Then, D^(x), 
^vXx) ^eO^) are the interpolating fields of the deuteron and leptonic pair, respectively. 

The effective Lagrangian Eq. (|2.3|) defines the effective vertex of the low-energy nuclear 
transition p + p— >D + e + + z/ e 



z.M(p + p^D + e+ + z/ e ) = G Y g A M N G n ^^e;(k D )[u(k Uc )Y(l-r)v(k e+ )} 

x [u c (p 2 )^u( Pl )}, (2.4) 



where e*(/co) is a 4-vector of a polarization of the deuteron, u(k Uc ), v(k c +), u(j>2) and 
u(j>i) are the Dirac bispinors of neutrino, positron and two protons, respectively. 

In order to evaluate the contribution of low-energy elastic pp scattering we have to 
determine the effective vertex of the p + p — > p + p transition accounting for the Coulomb 
repulsion between the protons. For this aim we suggest to use the effective local four- 
nucleon interaction Eq. (|2.1| ) and take into account the Coulomb repulsion in terms of the 
explicit Coulomb wave function of the protons. This yields 

V pp ^ p (k',k) = G„ NN r pp (k') [«(p 2 )tWi)] M P2 h 5 u(Pi)]^ PP (k), (2.5) 

where ip pp {k) and Vppt^') are the explicit Coulomb wave functions of the relative move- 
ment of the protons taken at zero relative radius, and k and k! are relative 3-momenta 
of the protons k = (pi — p*2)/2 and k' = (p{ — p 2 ')/2 i n the initial and final states. The 
explicit form of ip pp {k) we take following Kong and Ravndal [8] (see also [19]) 

^ PP (A;) = e-/ 4 ^r(l + ^-), (2.6) 

where re = l/M^a = 28.82 fm and a = 1/137 are the Bohr radius of a proton and the 
fine structure constant. The squared value of the modulo of ip pp {k) is given by 

| W *)|» = ^ W = ] i_ 7JE L_, (2.7) 

where Co{k) is the Gamow penetration factor [1,2,19]. We would like to emphasize that 
the wave function Eq.( f2~lf ) is defined only by a regular solution of the Schrodinger equation 
for the pure Coulomb potential [19]. 



By taking into account the contribution of the Coulomb wave function and summing 
up an infinite series of one-proton loop diagrams the amplitude of the solar proton burning 
can be written in the form 



iM(p + p^ D + e + + z/ e 



G v g A M N GWn fe e* u (k B ) [u(k u J^(l - 7 5 MM] F 
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(2. 



Mk-P-P-Q 1 m n -p-q 



where P = Pi + P2 = (2 y k 2 + M^, ) is the 4-momentum of the pp-pair in the center 
of mass frame; Q = aP + bK = a {jp\ + p?) + b (p\ — p?) is an arbitrary shift of virtual 
momentum with arbitrary parameters a and b, and in the center of mass frame K = 
Pi — Pi — (0,2 k) [14]. The parameters a and b can be functions of k. The factor J 7 ^ 
describes the overlap of the Coulomb and strong interactions [10]. It is analogous the 
overlap integral in the PMA [5]. We calculate this factor below. 

The evaluation of the momentum integral runs the way expounded in [14] . Keep only 
leading contributions in the large Nq expansion [13,14] we obtain 



J^\^p(\p+Q\)\ 2 ^ \r 



i 



i 



M^-p-P-Q 1 M N -p-Q 



-8 a (a + 1) Ml + 8 (b 2 -a(a + l))k 2 -i 8vrM N k \ip pp {k)\ 2 
-8 a (a + 1) Ml + 8 (b 2 -a(a + l))k 2 -i 8vr M N k C${k). 



(2.9) 



Substituting Eq.([0| in Eq.(gl| we get 
0/f(p + p ^D + e + + u e ] 



GV 9aMj<j Gttnn — - •?%! 



4vr 2 pp 
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2kr c 



l-a(a+l)^M 2 + ^(b 2 -a(a + l))k-- -, 
2n l 2ix l 



nn-^n 



2ty 



kC 2 (k) 



(2.10) 



In order to reconcile the contribution of low-energy elastic pp scattering with low-energy 
nuclear phenomenology [19] we should make a few changes. For this aim we should rewrite 
Eq. (p.lO|) in more convenient form 
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(2.11) 



(2.12) 



where a (k) is a pure Coulomb phase shift. 

Now, let us rewrite the denominator of the amplitude Eq.( |2.11|) in the equivalent form 
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re 



1 - i a PP r ; P fc2 + 7^ ^fcrc) + ^p P fc ^(^) 



(2.13) 



where we have denoted 

1 
Z 



1 9 PP PP 



+ cos ctq (A;) 



1 — a(o + 1)- 



h(2kr c ) 



G^nnMi 
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27T 2 
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(2.14) 



Here Z is a constant which will be removed the renormalization of the wave functions 
of the protons, a e pp = (-7.8196 ± 0.0026) fm and r e pp = 2.790 ± 0.014 fm [20] are the S- 
wave scattering length and the effective range of pp scattering in the 1 So-state with the 
Coulomb repulsion, and h(2krc) is defined by [19] 



h(2krc) = -7 + en(2kra) + £ „ (1+teW) - 



(2.15) 



The validity of the relations Eq. ( p.!4| ) assumes the dependence of parameters a and b on 
the relative momentum k. 

After the changes Eq. (|2. 1 lp -Eq. (|2.14| ) the amplitude Eq. (|2.10|) takes the form 



iM(p + P -> D + e + + u e ) = G y g A M N G ffNN ^ e*(k D ) [u(kj^(l - ^ 5 )v(k e+ )\ T. 
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pp 



X 
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Z [« c (p 2 ) 7 5 «(pi)]. 



(2.16) 



>'c 



Following [14] and renormalizing the wave functions of the protons \fZu{jpo) - > u(p2) and 
^u(pi) —> u(pi) we obtain the amplitude of the solar proton burning 



3g 



v 



iM(p + p - D + e + + u e ) = G v g A M N G nNN -^ e* (k D ) [u(K c )Y(l ~ 7>(M] T. 
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(2.17) 



re 



where we have introduced too an universal form factor [14] 

^ 2) = TTTfiF (2 ' 18) 

describing a spatial smearing of the deuteron coupled to the NN system in the 1 So-state at 
low energies; r D = I/V^d^n = 4.315 fm is the radius of the deuteron and e^ = 2.225 MeV 
is the binding energy of the deuteron. 

The real part of the denominator of the amplitude Eq. (|2.17| ) is in complete agreement 
with a phenomenological relation [19] 



ctgSlJk) 
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(2.19) 



describing the phase shift 5p P {k) of low-energy elastic pp scattering in terms of the S- 
wave scattering length a pp and the effective range rp* . As has been pointed out [19] the 
expansion Eq.( |2.19|) is valid up to T pp < 10 MeV, where T pp = k 2 /Mn is a kinetic energy 
of the relative movement of the protons. 

Thus, we argue that the contribution of low-energy elastic pp scattering to the am- 
plitude of the solar proton burning is described in agreement with low-energy nuclear 
phenomenology in terms of the S-wave scattering length a e and the effective range r p 
taken from the experimental data [20]. This takes away the problem pointed out by Bah- 
call and Kamionkowski [17] that in the RFMD with the local four-nucleon interaction 
given by Eq.( [2.1|) one cannot describe low-energy elastic pp scattering with the Coulomb 
repulsion in agreement with low-energy nuclear phenomenology. 

Now let us proceed to the evaluation of T p . For this aim we should write down the 
matrix element of the transition p + p^D + e + + i/ e with the Coulomb repulsion. The 
required matrix element has been derived in Refs. [11,14] and reads 

iM c {v + P ->• D + e+ + v e ) = 

= G Y g A M N G nNN ^ C {k) [u(k v ,)-f{l - 7 5 MM] <(&d) 

x {-[u^p 2 ) lal \( Pl )] J^(k D , h) - [u c { P 2)l h u{Pi)\ Jc U (kB, kt)}, (2.20) 

where kn and ke are 4-momenta of the deuteron and the leptonic pair, respectively. The 
structure functions J^ afil/ (k-D,ke) and J^ik^^kp) are determined by [11,14] 
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where q — p+ [kt — k^)j2. 



For the subsequent analysis it is convenient to represent the structure functions in the 
form of two terms 
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The decomposition is caused by the change 
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(2.23) 



where the first term gives the contribution to the SS part of the structure functions 
defined by strong interactions only, while the second one vanishes at re — > oo ( or a — > 0) 
and describes the contribution to the SC part of the structure functions caused by both 
strong and Coulomb interactions. 

The procedure of the evaluation of the structure functions Eq.( |2.21|) and Eq. (|2.22|) 
has been described in detail in Ref. [11,14]. Following this procedure we obtain J 7 ^ in the 
form 
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The integral can be estimated perturbatively The result reads 
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The numerical value of iJ-Lr is 
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(2.26) 



The contribution of the Coulomb field Eq. fl2.26p inside the one-nucleon loop diagrams is 
found small. This is because of the integrals are concentrated around virtual momenta 
of order of Mn which is of order Mn ~ Nc in the large Nc expansion [12]. For the 
calculation of the astrophysical factor 5 , pp (0) we can set J 7 ® = 1. 



3 Astrophysical factor for solar proton burning 



The amplitude Eq. fl2.17l ) squared, averaged over polarizations of protons and summed 
over polarizations of final particles reads 
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where m e = 0.511 MeV is the positron mass, and we have used the relation g^ = 
2n 2 Q D Ml 

In the low-energy limit the computation of the traces yields 
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where we have neglected the relative kinetic energy of the protons with respect to the 
mass of the proton. 

Substituting Eq. Q in Eq. Q we get 
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The integration over the phase volume of the final De + z/ c -state we perform in the non- 
relativistic limit 



d 3 k D d 3 k e + d 3 k Vc 

(2rr) 3 2E D (2n) 3 2E e+ (2n) 3 2E v 



(2n) A 5^(k D + k e -pi-p 2 ) \E c+ E Uc - h e+ ■ fc„. j 



1 



32vr 3 M, 



W+T D 



N J m e 



(W + T p 



E 2 + -m 2 E e+ (W + T pp -E e+ ) 2 dE c+ = ^ QQ J™ J /(£), (3. J) 



where W = £ D - (M n - M p ) = (2.225 - 1.293) MeV = 0.932 MeV and f = m e /(W + T pp ) 
The function /(£) is defined by the integral 
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and normalized to unity at £ = 0. 

Thus, the cross section for the solar proton burning is given by 
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MxT pp C*(JM N T pp ) 



Spp(T pp ) -7r/rcy M N T] 



pp 



T 



(3.6) 



pp 



The astrophysical factor S pp (T pp ) reads 
9g 2 A G 2 v Q D M? 



5 P p(Tj 



pp ' 



Q ■ 



1280tt 4 



G 2 NN (IU + T pp ) 5 /| 



m r 



W + T X 



PP/ 



^(MnTj 



pp; 



1 



1 - - a e r; p M N T pp + -E h(2r c JM N T, 



pp. 



+ o 



y^NTppCjCJMNTpp) 



(3.7) 



2 "" *" r " re 

At zero kinetic energy of the relative movement of the protons T pp = the astrophysical 



factor S'pp(O) is given by 
^pp(O) 



a 



dglGlQvM* q2 



1280tt 4 



- ff '/f 



4.08 x 10 _25 MeVb. 



(3.1 



The value Sp P (0) = 4.08 x 10~ 25 MeVb agrees good with the recommended value 5 PP (0) = 
4.00 x 10~ 25 MeVb [4]. Insignificant disagreement with the result obtained in Ref.[ll] 
where we have found <S , PP (0) = 4.02 x 10~ 25 MeVb is due to the new value of the constant 
g A = 1.260 -> 1.267 [18] (see Ref.[13]). 

Unlike the astrophysical factor obtained by Kamionkowski and Bahcall [5] the astro- 
physical factor given by Eq. ( |3.8j ) does not depend explicitly on the S-wave scattering wave 
of pp scattering. This is due to the normalization of the wave function of the relative 
movement of two protons. After the summation of an infinite series and by using the 
relation Eq. ( |2.19| ) we obtain the wave function of two protons in the form 

j5; v (k) sms; p (k) 

-a; p kC (kY 

that corresponds the normalization of the wave function of the relative movement of two 
protons used by Schiavilla et al. [6]. For the more detailed discussion of this problem we 
relegate readers to the paper by Schiavilla et al. [6] 1 . 

Unfortunately, the value of the astrophysical factor S , pp (0) = 4.08 x 10 -25 MeVb does 
not confirm the enhancement by a factor of 1.4 obtained in the modified version of the 
RFMDinRef.[14]. 



^pp(*0 



(3.9) 



4 Neutrino disintegration of the deuteron induced by 
charged weak current 

The evaluation of the amplitude of the process f e + D^e~+p + p has been given in 
details in Ref. [10]. The result can be written in the following form 

Gy 3(?v 



iM(p + p -»■ D + e + + v e ) = g A M h 



V2 2n< 



G.NN e!(M [u(M 7 ^l - 7 >(^ e )] ^' 



c 
ppe- 



■"■See the last paragraph of Sect. 3 and the first paragraph of Sect. 5 of Ref. [6]. 
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x J^ [<?(p 2 ) 7 5 u(pi)] Mk 2 ), (4.1) 

1 - g a l P r lp k2 + ^h(2kr c )+ i a c pp kC 2 Q (k) 

where J 7 * is the overlap factor which we evaluate below, and i<b(fc 2 ) is the universal 
form factor Eq. fl2.17|) describing a spatial smearing of the deuteron [14]. 

The amplitude Eq.( |4.1|) squared, averaged over polarizations of the deuteron and 
summed over polarizations of the final particles reads 

2 ,^6 144G v (5d n2 n-e 12 7T.2/r2> 



|.M(i/. + D-e-+p + p)| 2 = a lJW5 f^Gi^\F°\ 2 F^)F(Z,E, 



71 



x ; ,c C ° ik) (e c -E Vc - \k c - ■ k Ue ), (4.2) 



1 " 2°*^ + i^ ^ 2 M + K P Y k2 C 4 (k) 



3 



where F(Z, E e - is the Fermi function [21] describing the Coulomb interaction of the 
electron with the nuclear system having a charge Z. In the case of the reaction z/ e + D 
— > e + + p + p we have Z = 2. At a 2 Z 2 <C 1 the Fermi function F(Z, E e -) reads [21] 

F ^ E --)= r^=W- (4 ' 3) 



where rj e - = Za/v e - = ZaE e -/JE 2 _ — m 2 _ and v e - is a velocity of the electron. 

The r.h.s. of Eq.( |4.2| ) can be expressed in terms of the astrophysical factor 5 , pp (0) for 
the solar proton burning brought up to the form 
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o12k_2 r „M 3 \T e 

z 07r 'C lvl ti Kppe-I P 2^2> 

: ' i>iW; Vtr, + m 5 \T e I 2 



|^ (z/c + D ^e- + p + p)| 2 = 5 pp (0)— — ^—--pEl-F^^F^,^-) 



p P i 



x = ^^ (e c -E u - -k c - • £ V (4.4) 

"l " l^^k 2 + — h (2krc)] 2 + WtfCi(k) V 3 ' 

We have used here the expression for the astrophysical factor 



PPI 12807T 4 r c ^NNKppI "V 



^pp(O) = Toll N Gl NN \F;J 2 mln De+UB , (4.5) 



where m e = 0.511 MeV is the electron mass, and Qr>e+u c = (W/m e ) 5 f(m e /W) = 4.481 at 
W = 0.932 MeV. The function f(m e /W) is defined by Eq.Q. 

In the rest frame of the deuteron the cross section for the process u e + D — ► e~ + p 
+ p is defined as 

< D ^ = W. / W". + D-e- + p + p)|» 

i (2.)* «W( te + *;, - Pl - P2 - k,-) ^g^ p^ p^,^ . («) 

where E 1 ^, £i, £2 and _E e - are the energies of the neutrino, the protons and the electron. 
The abbreviation (cc) means the charged current. The integration over the phase volume 
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of the (ppe 
deuteron, 

1 
2 



-state we perform in the non-relativistic limit and in the rest frame of the 



d 3 pi d 3 p2 



d 3 k P 



(2n) 4 6 (4 \k D + K c -p 1 -p 2 -k e -) 



(2n) 3 2E 1 (2tt) 3 2E 2 (2tt) 3 2£ ( 

C 2 UM N T pp ) F 2 (M N T PP ) F{Z, E c - 



rc 



1 - 7ra e pp r; p M N T pp + -E£ h(2r c JM N T pp ) + (a; p ) 2 M N T pp C 4 ( V 'M N T pp ) 



E e -E Uc — —k e - ■ k Uc 



E Pb Mn (E th ^ /2 (2m^ /2 1 
128tt3 \M n ) \E th ) El 

T P . 



dT e - dT pp 5(E Uc - E th - T e - - T pp ) yT c - T pp I 1 + ^- 1 J 1 + g- 



C^(JM N T PP ) F 2 (M N T PP ) F(Z, £ e 



re 

xV2 / _ \3/J 



1 - 2 a PP r ;p M N T p P + ^ h(2r c ^/M N T pp ) + (a c pp ) 2 M N T pp C 4 ( V 'M N T I 
^M* [E th X -(2m. 



pp; 



128tt 3 \M N / \E, 



tb 



(j/-l) fippe-d/), 



(4.7) 



where T e - is the kinetic energy of the electron, E t h is the neutrino energy threshold of 
the reaction u e + D — >■ e~ + p + p, and is given by E t h = En + m e — (M n — M p ) = 
(2.225 + 0.511 - 1.293) MeV = 1.443 MeV. The function tt ppe -{y), where y = E u jE th , is 
defined as 



^ P pe- (y) = I dxy/x(l-x)(l + ^(y - 1)(1 - x)j Jl 



+ £<»-W-*) 



C^M N E th (y - 1) x) F 2 (M N £ th (y - 1) x) F(Z, m c + E th (y - 1) (1 - x)) 
1 " TTCC^N^th (y-l)x + -& h(2r c JM N E th (y-l)x) 



r c 



+ (< p ) 2 M N E th (y - 1) xC%(JM N E th (y-l)x)\, 



(U 



where we have changed the variable T pp = (E Uc — E th ) x. 
The cross section for u e + D — > e~ + p + p is defined 

7/2|jre |2 



3/2 



<f{E Va 



TT^Dc+i/o V E th 



\2m e 



p 12 
^ PPI 



3.69 x 10 5 S pp (0)J5^(i/- l) 2 n ppo -(y), 

K ppl 



(y-l^fippe-^) 



where S'pp(O) is measured in MeV cm 2 . For S pp (0) = 4.08 x 10 49 MeV cm 2 Eq.i 
cross section a^ D (E Uc ) reads 

<f(Ej = 1.50 l -^4- (V - !) 2 <W (v) ID" 43 cm 2 . 



I T e 1 2 
K ppl 



(4.9) 
) the 

(4.10) 
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In order to make numerical predictions for the cross section Eq. (|4.10| ) we should evaluate 
the overlap factor .7+™- ■ This evaluation can be carried out in analogy with the evaluation 
of J 7 * By using the results obtained in Ref.[10] we get 



t: 



ppe- 



1 + 



32 
9~ 



dvv z 



^—air/Av 



r 1- 



ia 
2v 



- 1 



The perturbative evaluation of the integral gives 

'13tt .137' 



PPe" 



1 — a 



54 



27 



(l +u 2)5/2 16(1 + ^)3/2 



+ Q(a' 



(4.11) 



(4.12) 



differs slightly from unity as well as the overlap factor T* 



Thus, the overlap factor J 7 ' 

of the solar proton burning. The ratio of the overlap factors is equal to 



pp 



I <pe 12 

l-^ppe-l 

I -re 12 
K PP I 



2tt 



a 



0(a 2 



l + (-1.53x 10 



-2\ 



1. 



(4.13) 



Setting l-^ppe-P/l-^ppI 2 = 1 we can make numerical predictions for the cross section 
Eq. (|4.10 ) and compare them with the PMA ones. 

The most recent PMA calculations the cross section for the reaction u e + D — > e~ + 
p + p have been obtained in Refs. [22,23] and tabulated for the neutrino energies ranging 
over the region from threshold up to 160 MeV. Since our result is restricted by the neutrino 
energies from threshold up to 10 MeV, we compute the cross section only for this energy 
region 



o. 



<f{E Ve 


= 4 MeV) 


= 2.46 (1.86/1.54) x 10" 43 cm 2 


<f{E Vc 


= 6 MeV) 


= 9.60 (5.89/6.13) x 10~ 43 cm 2 


<f(E Uc 


= 8 MeV) 


= 2.38 (1.38/1.44) x 10" 42 cm 2 


■ u c f(Eu c = 


= 10 MeV) 


= 4.07(2.55/2.66) x 10" 43 cm 2 



(4.14) 



where the data in parentheses are taken from Refs. [22] and [23], respectively. Thus, on 
the average our numerical values for the cross section a^ (E Ue ) by a factor of 1.5 are 
larger compared with the PMA ones. 

Our predictions for the cross section Eq. fl4.14|) differ from the predictions of Ref.[14]. 
This is related to (i) the value of the astrophysical factor which is by a factor 1.4 larger 
in Ref . [14] and (ii) the form factor describing a spatial smearing of the deuteron which is 
F 2 (fc 2 ) is this paper (see Ref. [13]) and F D (k 2 ) in Ref. [14]. 



5 Astrophysical factor for pep process 

In the RFMD the amplitude of the reaction p + e~ + p — > D + u e or the pep-process is 
related to the effective Lagrangian Eq.( ^.3[ ) and reads 



3g 



v 



iM (p + e" -b p - D + u e ) = GV ff A M N GW N ^ e* (A*>) [u(K c )^(l - i 5 )u(k e -)} F\ 



c 

pp 



Co(fc) 



5 

re 



- al p rl v k 2 + ^ h(2kr c ) +i a c pp k C 2 (k) 



[^(p 2 ) 7 5 n( Pl )]F D (fc 2 ), 



(5.1) 
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where we have described low-energy elastic pp scattering in analogy with the solar proton 
burning and the neutrino disintegration of the deuteron. 

The amplitude Eq. ([5.1|) squared, averaged and summed over polarizations of the in- 
teracting particles is defined 

|A*(p + e-+p-D + I / B )|2 = G\ giM« G 2 NN HZ^ \j^\* F *(k 2 ) F(Z, E c -) 

7T Z 

C2 ° {k) ( E e+ E Ve - l -k e+ -t\ (5.2) 



1 ^i^ 2 ^M2M] 2 H-K P W(fc) 1 ° 3 ' 

where F(Z, E e -) is the Fermi function given by Eq.( f4.3| ). 

At low energies the cross section cr pep (T pp ) for the pep-process can be determined as 
follows [24] 

j 1 r d ]c — -*■ f 

(2vr) 4 ^(A; D + ^ - Pl - p 2 - fc e -) , ^ D , ^ , (5.3) 

where g = 2 is the number of the electron spin states and v is a relative velocity of the 
protons. The electron distribution function n(k e -) can be taken in the form [21] 

n(k e -) = e jy - T --/ kTc , (5.4) 

where k = 8.617 x 1CT 11 MeVK^ 1 , T c is a temperature of the core of the Sun. The 
distribution function n(k e -) is normalized by the condition 

9 J 7^r~p n (^-) =ne "' ( 5 - 5 ) 

where n e - is the electron number density. From the normalization condition Eq.( |5.5| ) we 
derive 

e* = 4n3ne \ /2 . (5.6) 

(2vrm e fcT c ) 3/2 V ; 

The astrophysical factor S , pep (0) is then defined by 
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2ir flj) e + 

For the ratio S pep (0) / S pp (0) we obtain 

S pep (0) _ 2 3 /y/ 2 ( aZn e - \ (E th V [m[J 7 [^\ , . 

S PP (0) /pp(0) ^ ml ) \mj \jkT c \ V kT c )' l ' } 

We have set / pp (0) = VLd c + Uc /?>0 = 0.149 [21] and the function I(x) having been introduced 
by Bahcall and May [21] reads 

/(») = / *""" . r . (5.9) 

o 1 — e ' v 
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— Af-V^ fd%-e- T <-/^F(Z,E, 

c+v c m%\m e J J 



S pep (0) = 5 pp (0) ^ — ( — I * V ld%-e-^-/ k ^F(Z,E e -). (5.7) 



The relation between the astrophysical factors S , pep (0) and 5 , pp (0) given by Eq. (|5.8| ) is in 
complete agreement with that obtained by Bahcall and May [21]. The ratio Eq.( |5.8| ) does 
not depend on whether the astrophysical factor 5 , pp (0) is enhanced with respect to the 
recommended value or not. 



6 Conclusion 

We have shown that the contributions of low-energy elastic pp scattering in the x So-state 
with the Coulomb repulsion to the amplitudes of the reactions p + p— >D + e + + ^ e , 
v e + D — >• e~ + p + p and p + e~ + p^D + i/ e can be described in the RFMD in 
full agreement with low-energy nuclear phenomenology in terms of the S-wave scattering 
length and the effective range. The amplitude of low-energy elastic pp scattering has 
been obtained by summing up an infinite series of one-proton loop diagrams and the 
evaluation of the result of the summation in leading order in the large Nc expansion. 
This takes away fully the problem pointed out by Bahcall and Kamionkowski [17] that in 
the RFMD with the effective local four-nucleon interaction Eq. (|2.1|) one cannot describe 
low-energy elastic pp scattering in the l So-state with the Coulomb repulsion in agreement 
with low-energy nuclear phenomenology. 

The obtained numerical value of the astrophysical factor S'pp(O) = 4.08 x 10~ 25 MeVb 
agrees with the recommended value 5 , pp (0) = 4.00 x 10~ 25 MeVb and recent estimate 
S'pp(O) = 4.20 x 10~ 25 MeVb [9] obtained from the helioseismic data. 

Unfortunately, the value of the astrophysical factor 5 , pp (0) = 4.08 x 10~ 25 MeVb does 
not confirm the enhancement by a factor of 1.4 obtained in the modified version of the 
RFMD in Ref . [14] which is not well defined due to a violation of Lorentz invariance of the 
effective four-nucleon interaction describing N + N — > N + N transitions. This violation 
has turned out to be incompatible with a dominance of one-nucleon loop anomalies which 
are Lorentz covariant. 

The cross section for the neutrino disintegration of the deuteron has been evaluated 
with respect to 5 , pp (0). We have obtained an enhancement of the cross section by a 
factor of order 1.5 on the average for neutrino energies E Vc varying from threshold to 
E Uc < 10 MeV. It would be important to verify our results for the reaction v e + D — > e~ 
+ p + p in solar neutrino experiments planned by SNO. In fact, first, this should provide 
an experimental study of 5 , pp (0) and, second, the cross sections for the anti-neutrino 
disintegration of the deuteron caused by charged u e + D — »e + +n + n and neutral 
u e + D — >z/ e + n + p weak currents have been found in good agreement with recent 
experimental data obtained by the Reines's experimental group [26]. 

The evaluation of the astrophysical factor Sp^O) for the reaction p + e~ + p — > D 
+ u e or pep-process in the RFMD has shown that the ratio S pep (0) / S pp (0) , first, agrees 
fully with the result obtained by Bahcall and May [21] and, second, does not depend on 
whether 5 , pp (0) is enhanced with respect to the recommended value or not. 

Concluding the paper we would like to emphasize that our model, the RFMD, con- 
veys the idea of a dominant role of one-fermion loop (one-nucleon loop) anomalies from 
elementary particle physics to the nuclear one. This is a new approach to the description 
of low-energy nuclear forces in physics of finite nuclei. In spite of almost 30 year's history 
after the discovery of one-fermion loop anomalies and application of these anomalies to 
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the evaluation of effective Lagrangians of low-energy interactions of hadrons, in nuclear 
physics fermion-loop anomalies have not been applied to the analysis of low-energy nu- 
clear interactions and properties of nuclei. However, an important role of NN fluctuations 
for the correct description of low-energy properties of finite nuclei has been understood 
in Ref.[16]. Moreover, NN fluctuations have been described in terms of one-nucleon loop 
diagrams within quantum field theoretic approaches, but the contributions of one-nucleon 
loop anomalies have not been considered in the papers of Ref.[16]. 

The RFMD strives to fill this blank. Within the framework of the RFMD we aim 
to understand, in principle, the possibility of the description of strong low-energy nu- 
clear forces in terms of one-nucleon loop anomalies. Of course, our results should be 
quantitatively compared with the experimental data and other theoretical approaches. 
Nevertheless, at the present level of the development of our model one cannot demand at 
once to describe, for example, the astrophysical factor 5 , pp (0) with accuracy better than 
it has been carried out by Schiavilla et al. [6], where only corrections not greater than 
1% are allowed. It is not important for our approach at present. What is much more 
important is in the possibility to describe without free parameters in quantitative agree- 
ment with both the experimental data and other theoretical approaches all multitude of 
low-energy nuclear reactions of the deuteron coupled to nucleons and other particles. In 
Ref. [13] we have outlined the procedure of the evaluation of chiral meson-loop corrections 
in the RFMD. The absence of free parameters in the RFMD gives the possibility to value 
not only the role of these corrections but also the corrections of other kind mentioned 
recently by Vogel and Beacom [25]. 

The justification of the RFMD within QCD and large Nq expansion [12] implies that 
one-nucleon loop anomalies might be natural objects for the understanding of low-energy 
nuclear forces. The real accuracy of the approach should be found out for the process of 
the development. 
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